Conditions for spontaneous homogenization of the Universe 

Krzysztof Bolejko* 

■ Steward Observatory, University of Arizona, Tucson, AZ 85721, USA 

^ ! and Nicolaus Copernicus Astronomical Center, Bartycka 18, 00-716 Warszawa, Poland 

^' William R. Stoeger 

Vatican Observatory Research Group, University of Arizona, Tucson, AZ 85721, USA 

O : May 18, 2010 

^^, Abstract 

$J_( , The present-day Universe appears to be homogeneous on very large scales. Yet when the casual 

Olj, structure of the early Universe is considered, it becomes apparent that the early Universe must have 

been highly inhomogeneous. The current paradigm attempts to answer this problem by postulating 

the inflation mechanism However, inflation in order to start requires a homogeneous patch of at least 

^ ' the horizon size. This paper examines if dynamical processes of the early Universe could lead to 

^\ ' homogenization. In the past similar studies seem to imply that the set of initial conditions that 

c7^ . leads to homogenization is of measure zero. This essay proves contrary: a set of initial conditions for 

^Z^ ' spontaneous homogenization of cosmological models can form a set of non-zero measure. 
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The present-day Universe appears to be homogeneous on very large scales. The main argument comes 
from the isotropy of the cosmic microwave background radiation (CMB). The Ehlers-Geren-Sachs (EGS) 
theorem [T] and 'almost EGS theorem' |2] imply that if anisotropics in the CMB are small for all funda- 
mental observers then locally the Universe is almost spatially homogeneous and isotropic. In addition, 
if there were large inhomogeneities in the Universe they would manifest themselves in CMB tempera- 
ture fluctuations via the Ree-Sciama effect [31I31[S]. Moreover, the enormous success of the homogeneous 
Friedmann-Lemaitre-Robertson- Walker (FLRW) models in describing cosmological observations seems to 
provide further evidence for the large scale homogeneity of our Universe — nearly all sets of cosmological 
observations are successfully analyzed within the framework of homogeneous models, including supernova 
observations [51 [7], baryon acoustic oscillations [S] and the CMB [5]. 

On the other hand, when the casual structure of the early Universe is considered, it becomes apparent 
that the early Universe must have been highly inhomogeneous [101 E] ■ This raises a key question: how 
did the chaotic behavior expected in the very early Universe come to generate the Cosmos that on large 
scales looks so homogeneous? 

The current paradigm attempts to answer this question by postulating the inflation mechanism [121113] . 
Inflation not only appears to solve the horizon problem but also other problems of the standard cosmology 
like the origin of matter fluctuations (during inflation quantum fluctuations are enhanced and they become 
a source of primordial fluctuations) , magnetic monopole (while the particle physics predicts existence of 
magnetic monopoles or other exotic particles, these have not been observed as during inflation they were 
pushed outside the horizon), flatness problem (since the inflating region increases its size by factor ^ lO'^'^ 
the effective curvature decreases by ^ 10^^^), and many others. While alternative theorems such as cosmic 
stings or bounce cosmologies can deal with some of the problems they do not solve all the problems in 
such a graceful and simple way as inflation does |14| . 

However, after closer investigation it becomes apparent that inflation does not solve the problem of 
homogenization of the Universe. Rather, it explains how it is possible that regions which appear presently 
causally disconnected were casually connected at some very early stage of the evolution of the Universe. 
This is because the mechanism of inflation relies on explicit models that are homogeneous from the very 
beginning. The studies of inflation within inhomogeneous backgrounds proved that in order to start, 
inflation requires a homogeneous patch of at least the horizon size [151 HSl HZl [H] ■ 

Ironically, inflation which had originally been postulated to explain homogeneity, now relies on the 
assumption that the early Universe was homogeneous. This is what the lay people call a vicious circle. 
One may try to avoid the problem by pointing that there is always a probability that a homogeneous patch 
(even of a horizon size) existed in the Universe, or even by invoking the Anthropic Principle arguing that 
if inflation could not start then the galaxies and life could not form. However, the above arguments do not 
seem to be satisfactory — using similar arguments one could try to argue that the Universe was from the 
very beginning homogeneous and the initial perturbations were encrypted in it and then there would be 
no need for inflation at all. Given that we can study only one Universe, the problem of initial conditions 
of the Universe seems beyond the reach of modern science ^llj . Also, our current knowledge about matter 
and its properties under extreme conditions seems rather insufficient to study conditions of the early 
Universe which most probably should be described by the theory of quantum gravity. Nevertheless, the 
problem of homogenization should also be studied using classical methods. An elegant solution would 
be to prove some dynamical scenario (other than sheer luck) in which starting from arbitrary initial 
conditions the Universe would undergo spontaneous homogenization. On the other hand, if one can show 
that no dynamical process can homogenize the Universe before the start of inflation this would be a strong 
suggestion that initial conditions had to be tuned in order to achieve homogenization, or that there must 
be some other process that can explain current cosmological problems without inflation. 

The problem of dynamical homogenization of the Universe is not new. In 1968 Misner [19j proposed 
that the Universe started off in a chaotic state and that dissipation processes damped out inhomogeneities 
leading to a homogeneous state. Misner considered the Bianchi type I solution with viscosity. Among 
homogeneous and anisotropic models the Bianchi I models are the simplest and in the case of vanishing 
shear they reduce to the parabolic FLRW solution, and they are characterized by the following metric 



ds^ = -dt^ + X^{t)dx^ + Y^(t)dy^ + Z^{t)dz^. (1) 

It can be shown [20j that the function S'^{t) = XYZ evolves as 

where S is a constant. As seen if i ^ S, anisotropies decrease — consider a simple case oi Y = Z then 
the shear is equal to 

"'^3^ (lis)'- P) 

For t ^ T,, a'^ ^ a~^, where a{t) is the scale factor of the FLRW model, and S{t) > a{t). If instead of 

dust a viscous fluid is considered then at late times, as shown by Ellis |21) 

a^ ^ !^Piz^, (4) 

a° 

where A is the viscosity coefficient. This simple example shows that viscosity can induce an exponential 

decay of the shear and therefore be a source of isotropization. 

However, Collins & Hawking [5^ showed that among homogeneous and anisotropic models those 

that approach homogeneous and isotropic state are of measure zero. They also suggested that within 

inhomogeneous framework homogenization seems unlikely. On the other hand Bonnor and Tomimura 

P51 ^^ showed that within a wide class of initial conditions inhomogeneous models (e.g. Lemaitre- 

Tolman and Szekeres geometries) which evolve into the homogeneous FLRW models can exist. In other 

words, there exist inhomogeneous solutions with only decaying modes. As later confirmed by Silk [5S] 

and Plebahski & Krasiiiski [55] such models are characterized by being free of curvature perturbation. 

One may think that models with only decaying modes form a set of measure zero, since one of arbitrary 

functions defining such models must be of specific form. However, there are some cases when perturbations 

are described only by the decaying modes. For example let us consider a 'comformal-Newtonian' metric 

ds2 = a^ [(1 + $)dr/2 - (1 - 2*)%dxMd^] , (5) 

For the ultra- relativistic matter or radiation-dominated Universe {p/p = 1/3) the adiabatic perturbations 
evolves as [13] 
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where x = kri/\/3, fc is a wavenumber, and 77 is a conformal time {dt = adrj). Since $ ~ —25p/p then 
the adiabatic perturbations of the ultra-relativistic matter are always described by decaying modes. Since 
at high temperatures matter should be ultra-relativistic, this shows that the early Universe could provide 
suitable conditions for homogenization. 

So far we have considered only simple examples, now we consider whether in more general inhomoge- 
neous case the set of initial conditions that lead to homogenization can be of non-zero measure. Let us 
assume that there is a unique timelike vector field u"^ which can be associated with the average velocity 
of matter [21j , then the tensor 

h\ = S\ + u'^Ub, (7) 

projects tensors into the surface orthogonal to m°. Further the velocity field can be decomposed [21] 

Ua;b = (^ab + (yab + -::hab0 - UaUb, (8) 



into shear {(Tab), rotation {uiab), expansion (9) and acceleration {iia — Ua^bU^)- The homogeneous solutions 
form a set of vanishing rotation, shear and acceleration. Therefore, a potential process of homogenization 
of the universe must result in reducing the amplitude of uJab-, cFab-, and Ua to some negligible values. Let 
us focus on the shear evolution which is given by |21] 
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hJhb^&fg~hJhb^U(f.g)-UaUb+i^aUJb + aJafb + -0<Tab + hab{--:;i^'^--;:O-^ + :^u''c)--^^ = 0. (9) 

Since the propagation equation for cjab and Ua are of similar complexity it is hard to say, without detailed 
numerical calculations how the Universe can homogenize itself. However, if the Universe is expansion- 
dominated, i.e the term 9aab dominates in the above equation, then the shear decays. In the early Universe 
the expansion rate is indeed of high amplitude, but so the other terms can be. To build up our intuition 
about dominance of Oaab let us consider a few simple examples. For simplicity we will consider space-times 
that are locally rotationally symmetric (LRS) . In such a case at each point there exists a unique preferred 
spatial direction which constitutes a local axis of symmetry [27]. Therefore, there exists a spacelike vector 
field such that 

eV, = 0, e,e^ = 1 ^ e'ei,^ = 0. (10) 

Using this vector field one can define a trace-free symmetric tensor 

Due to the LRS property all covariantly defined spacelike trace- free symmetric tensors which are orthogonal 
to u" have to be proportional to Cij. Therefore, the shear and the electric part of the Weyl tensors can 
be written as follows 
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where a = ^a^ja^i and E — ^E^jE^i. 

Let us now consider the Lemaitre-Tolman model {p = = lo ~ u). The metric of the Lemaitre-Tolman 
model is [28, 29J 

R ^ 

ds^ = -dt^ + I^TT^dr^ + R^{t, r)(cW2 + sin^ Mip^), (13) 

where E{r) is an arbitrary function, i?,^ — dR/dr, and R{t,r) obeys 

i?„2 = 2£; + ^ + |i?^ (14) 

where R,t = R_aU°', M{r) is another arbitrary function and A is the cosmological constant. Below it will 
be assumed that A = 0, first because the cosmological constant is negligible in the early universe, and 
second because we are interested in decelerating cosmology (in contrast to accelerated phase of inflation). 
The interpretation of the functions M{r) and E{r) is as follows: M is the gravitational mass contained 
within the comoving spherical shell at any given r, while E is the energy per unit mass of the particles on 
that shell. Moreover, E determines the space curvature at each r-value. For the Lemaitre-Tolman model 
eq. ([9]) reduces to 



where p = |^ and p = ^ (in the FLRW hmit p = p). Since 9 = 2T + P and a = ^{P - T) — where 

P = R' /R' is the parallel Hubble parameter, and T = R/ R is the tangential Hubble parameter — the 
second term on the rhs dominates over the first one. On the other hand the last term is 

It is therefore easy to notice that the parabolic Lemaitre-Tolman model {E = 0) undergoes homogeniza- 
tion. The same is true if the second and third terms of the rhs are much smaller than 2T{P — T). Let us 
consider, however, a stronger condition 

The above condition has an intuitive simple, geometric interpretation: it is equivalent to '^TZ <^ TZ (where 
■^7?. is the Ricci scalar of 3-space, and TZ is the Ricci scalar of space-time). Therefore, if only p7)) holds 
universe homogenize itself. 

Surprisingly, eqs. p3|) and p7|) are literally the same for the quasi-spherical Szekeres model (30] 
(which do not have any Killing vectors) and also for its generalization in the Roy-Singh model |31[ [32] 
(where bulk and shear viscosity are of such a special form so that u°- = 0). For the Roy-Singh model 
Eab — ■^{p ~ p)eafc + ^T^ab, and thus the viscosity cancels out in (j9]). This shows that the spherically 
symmetric dust solutions of the Einstein equations are not the only one that undergo homogenization 
under condition pT]). 

Keeping this in mind let us consider a more general case, when the source of the gravitational field is a 
fluid with inhomogeneous pressure, viscosity, and acceleration. Let us however keep spherically symmetry 
and iOab = 0. The most general form of spherical symmetric metric is 

ds^ = -e^di^ -I- e-^dr^ -|- R^M^ + R^ sin^ ddip^ . (18) 

In this case eq. © reduces to 

where V = ^e^^ Ua" - \A'B' -(- \A'^ - \^A'\ . Given that A'e-^ = {2E)/{RR') the following con- 
ditions 

l^«^' ^«^' |^«^' ('«) 

are sufficient to ensure that (2/3)9a term is larger than others. Although under these conditions shear 
vanishes, this is not a sufficient condition for homogenization, as there are solutions of the Einstein 
equations that are shear-free but inhomogeneous (see [33l[34]). However, if in addition it is assumed that 
the rate of change of entropy is positive then, as shown by Ellis [H] , this implies 

TTob = — AcTob, (21) 

where A > is a viscosity coefficient. So under the above condition vanishing shear implies vanishing 
viscosity and hence ^U^ leads to homogenization. 

Let us consider what these relations (pn|) imply. They imply that the function E is confined or small 
in comparison to contribution from density. However, such a geometric explanation does not provide 
physical insight. In the case of vanishing heat flux, the conservation equations {T'^^-i, = 0) reduce to 

^/ -p' + -^{Xcr)' + 2V3X<jR'/R 

— = ^ . (22) 
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Combining this with expression for T) and demanding that all terms remain smaller than p implies 

which shows that although the gradients of pressure, shear and viscosity can tend to infinity when i — > 
they cannot be unbound and must be smaller than p^'^. Given the fact that when matter is highly dense, 
its equation of state becomes degenerate (being weakly depend on temperature) it might happen that 
pressure and its gradients do not increase so rapidly as density as t — > 0. If this is indeed the case then 
conditions (j20p may be more probable than it first appears. 

Summarizing, it seems clear there are conditions in the early Universe that can lead to spontaneous 
homogenization and which are not of measure zero. This is in contrast to study of decaying modes in 
the Lemaitre-Tolman or FLRW models. The difference is that it is not required that the model becomes 
homogeneous aat—i'oo, but rather in the early stages of its evolution. If the density is large and gradients 
of pressure and viscosity not extreme, then relations pO|) ensure homogenization. Thus as long as the 
density is of very large amplitude the homogenization proceeds. It is interesting to notice that viscosity 
helps in keeping density at higher values. If heat flow is zero then the conservation equation T'^^-i, = 
reduces to 

p = -9{p + p) + Xa'^. (24) 

This relation shows that in the presence of viscosity (A<t^) density decreases more slowly than in the case 
of zero-viscosity. In addition as shown in Ref. |35| if gradients of pressure are present then the increase of 
density contrast in density waves (which can develop a shock when p = 0) is slower than in p' = case, 
and eventually the gradients lead to acoustic oscillations. Intuitively, one can imagine that if viscosity 
is present then it is harder to induce such oscillations thus keeping the density contrast at lower values. 
This all leaves an impression that dynamical processes of the early universe could lead to homogenization. 
Whether or not this took place in our Universe remains to be examined by detailed numerical calculations. 
Even if our Universe was not homogenized by some dynamical process, the main conclusion of this essays 
ramins unchanged: a set of initial conditions for spontaneous homogenization of cosmological models can 
form a set of non-zero measure. 
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